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Abstract 

We study the stability of holonomic quantum computations with respect to errors 
in assignment of control parameters. The general expression for fidelity is obtaned. 
In the small errors limit the simple formulae for the fidelity decrease rate is derived. 
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Holonomic approach to quantum computations was primarily proposed in the 
Ref. [1]. Optical holonomic quantum computer was proposed to realize this 
idea in a non-linear Kerr medium with the degenerate states of laser beams to 
be interpreted as qubits and logical gates to be realized by employing the exist- 
ing devices of quantum optics [2] . Another implementation of holonomic quan- 
tum computer based on trapped ions in optical cavity was discussed recently 
from the viewpoint of its resilience to control parameter errors (undesirable 
fluctuations of the laser beams amplitudes) [3]. Fidelity primarily proposed in 
the Ref. [4] is widely used as a measure of the stability of quantum systems. 
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Fidelity close to unity means stability of quantum system with respect to small 
perturbations of the Hamiltonian (or more generally evolution operator) [5]. 
Particulary it defines the resilience of quantum computations [6] . Another ap- 
proaches to determine the stability of quantum systems were also proposed 
(for instance see the Ref. [7]). It was demonstrated in particular for cavity 
model of holonomic quantum computer that this realization of logical gates 
provides fidelity and a success rate close to unity [8]. In this Letter we derive 
the general expression for the fidelity of any holonomic quantum computation. 
Taylor expansion for fidelity as function of the control errors magnitude for 
the general case of holonomic quantum computer is obtained in the most rel- 
evant limit of small errors. For this purpose the non-abelian Stokes theorem 
primarily utilized in connection with the confinement problem of quarks and 
gluons [9] is used. The simple formulae for the fidelity decrease rate when con- 
trol errors grow is extracted. It is worth to note that fidelity in particular for 
holonomic quantum computer is mathematically similar to the Wilson loop 
defined in the theory of strong interactions of elementary particles (quantum 
chromodynamics), see the Ref. [10]. Namely, the "area law" for Wilson loop 
means confinement of quarks and gluons, as well the similar behavior of fi- 
delity demonstrates the instability of quantum system. This corresponds to 
the connection between quantum chaos and confinement in the sense of the 
Ref. [7]. 

In holonomic quantum computer non-abelian geometric phases (holonomies) 
are used for implementation of unitary transformations (quantum gates) in the 
subspace C N spanned on eigenvectors corresponding the degenerate eigenvalue 
of parametric isospectral family of Hamiltonians 
F = {H(X) = U(\)H U(\) + }\ e M [!]• The A's are the control parameters and 
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M represents the space of the control parameters. The subspace C N is called 
quantum code (N is the dimension of the degenerate computational subspace). 
Quantum gates are realized when control parameters are adiabatically driven 
along the loops in the control manifold M. The unitary operator mapping the 
initial state vector into the final one has the form 0^ e^T^A 1 ), where / 
enumerates the energy levels of the system, (pi is the dynamical phase, R is 
the number of different energy levels of the system under consideration and 
the holonomy associated with the loop 7 G M is given by [11]: 



r 7 (A M ) =Pexp* J ApdXf,,. (1) 

7 

Here P denotes the path ordering operator, A^ is the matrix valued adiabatic 
connection given by the expression: 

(AOmn = / d 3 x^* m (x, A) A), (2) 



where integration goes over the spatial coordinates x and vp n (x, X), n = 1, N 
are basis functions of the corresponding eigenspace C N . Dynamical phase will 
be omitted bellow due to the suitable choice of the zero energy level [1]. We 
shall consider the single subspace (no energy level crossings are assumed). 

Fidelity for holonomic quantum computer in the general case can be written 
as follows: 



f = tr (pr;, 1 ^) = tr l^pP exp {1 j A^dX^ 



(3) 



where 70 is the adiabatic loop implementing the desirable quantum gate, 7' is 
the actual loop with some deviations from 70 due to practically unavoidable 
errors in assignment of control parameters, p denotes the density matrix of 
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the quantum state which stability is studied. The loop £7 is defined by the 
relation £7 = 7 ,_1 • 70. The -operation means that to obtain the contour £7 
we primarily go over the loop 70 in the straight direction and after that travel 
along the loop 7' in the opposite direction. For the definition and properties 
of this operation see Ref. [1]. The contour £7 is assumed to be closed i.e. to 
form a loop. 

Using the non-abelian Stokes theorem [9] is it easy to obtain the following 
expression for the fidelity (3): 



Here F xe = d x A Q — d e A x — i[A x , A e }_ is the curvature tensor in the space of 
the control parameters, d x = d/d\ x , [, ]_ denotes the commutator, SS is the 
surface spanned on the loop £7 and da xs is the projection of the infinitesimally 
small surface element of SS on the coordinate plane XQ of the control parameter 
space. In the formulae (4) and everywhere bellow the summation over the 
indices of the curvature tensor F xe is performed under the condition x > Q- 
Current point z is located on the loop £7 and point xq is the arbitrary one. In 
the Ref. [9] it was demonstrated that when the non-abelian Stokes theorem 
is applied the result does not depend on the particular form of SS and the 
point Xq is arbitrary if there are no monopole-like and string-like topological 
structures. This is the case assumed here. The general case deserves further 
investigation. 

The general expression (4) is not convenient for the practical use. To sim- 
plify it we consider more relevant limit of small errors. Namely, we assume 
that the errors in the assignment of control parameters <5A M are small and sat- 



/ = trpP exp % 



ss 



I 




(4) 
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isfy the restrictions \5\^\ < \5\ x 5\ p \ < (x > P)- The 

connection and curvature tensor are calculated in some point Ao defining 
the position of the small loop £7. The norm is defined as follows \\B\\ = 
sup{^J< ^\B+B\tp > I \ip >G C N , < ip\i/j >= 1}. Here B denotes both the 
connection matrix components and the curvature tensor matrix compo- 
nents F xp (x > p)- Then the following Taylor expansion in the point Ao G M 
can be obtained for the fidelity (4): 

/ = 1 + itrpF xg 5\ x 5\ e - trp[A p , F xg }-5\^5\ x 5\ e + (itrpA^F xe A u - 

- l -trp[F xs , A,A U ] + - l -trpF xe F^ 5\ x 5\ e 5\^8\„ + 0(5X 5 ). (5) 

Here the point A defines the position of the small loop £7 in the space of 
parameters M and [, ]+ denotes the anti-commutator. Under the stated re- 
striction on the control parameter errors higher order terms in Eq.(5) give less 
significant contribution to the fidelity then the terms having the lower order 
on 5\. It is seen that the Taylor expansion for fidelity does not have the linear 
on 5\ term. The cancellation of the first order terms was noticed for the par- 
ticular implementation of holonomic quantum computer on trapped ions [3]. 
We have demonstrated that this conclusion does not depend on particular 
realization of holonomic quantum computation. The reason is that the differ- 
ence between two vector parallel transported along two different infinitesimally 
small paths from one point to another is proportional to the area enclosed by 
these contours and it does not depend on their lengths. 

Equality of the fidelity to the unity means that the computations are stable. 
Namely, the result of computations is not changed, when the quantum gate is 
perturbed. Expression (5) shows that fidelity slightly deviates from unity for 
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small errors in assignment of control parameters. For zero errors fidelity equals 
unity as it should be. The rate of fidelity deviation from unity extracted from 
the second item of the expansion (5) is determined by the curvature tensor: 

= \tr P F^\ (6) 

where 5S^ U = 5\^5\ u is the "error area" in the plane \iv. We see that the 
less the curvature tensor leads the less the deviation of the fidelity from unity 
and therefore calculations are more stable. This result agrees with one of 
the Ref. [3], where for the particular realization of holonomic quantum com- 
puter on trapped ions it was demonstrated the one-qubit holonomic gate to 
be resilient to control errors for large values of the squeezing parameter (expo- 
nentially small values of the curvature tensor components). From the physical 
point of view this result can be understood as follows. The curvature tensor 
(as well as the area enclosed) determines the difference between two vector 
parallel transported along two different infinitesimally small paths from one 
point to another. Therefore the less the curvature means the less the differ- 
ence between state vectors parallel transported along loops 70 and 7'. Thus 
the computation is more stable. The condition for robust holonomic quantum 
computation is F^ v = for points belonging to the loop 7 . 

In conclusion, exploring the similarity between the mathematical apparatus 
used in non-abelian gauge field theories and one used for description of holo- 
nomic quantum computation we obtained the general expression for fidelity 
of holonomic quantum computer. Its Taylor expansion in the small control 
errors limit is derived and the simple formulae for the rate of the fidelity de- 
viation from the unity is obtained. Our general results are in agreement with 
ones obtained for the particular realization of holonomic quantum computer 
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on trapped ions. The formal mathematical analogy between the fidelity and 
Wilson loop was noted. The question about the relationship between the con- 
finement of the strongly interacting quarks and gluons and the stability of the 
corresponding systems is under the investigation. 
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